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ABSTRACT 

We consider the scenario of a magnetic field orthogonal to a front separating 
two media of different temperatures and densities, such as cold and warm neutral 
interstellar gas, in a 2-D plane-parallel geometry. A linear stability analysis is 
performed to assess the behavior of both evaporation and condensation fronts 
when subject to incompressible, corrugational perturbations with wavelengths 
larger than the thickness of the front. 

We discuss the behavior of fronts in both super-Alfvenic and sub-Alfvenic 
flows. Since the propagation speed of fronts is slow in the ISM, it is the sub- 
Alfvenic regime that is relevant, and magnetic fields are a significant influence on 
front dynamics . In this case we find that evaporation fronts, which are unstable 
in the hydrodynamic regime, are stabilized. Condensation fronts are unstable, 
but for parameters typical of the neutral ISM the growth rates are so slow that 
steady state fronts are effectively stable. However, the instability may become 
important if condensation proceeds at a sufficiently fast rate. 

This paper is the first in a series exploring the linear and nonlinear effects of 
magnetic field strength and orientation on the corrugational instability, with the 
ultimate goal of addressing outstanding questions about small-scale ISM struc- 
ture. 

Subject headings: instabilities — ISM: structure — MHD 

1. Introduction 

The interstellar medium (ISM) exhibits structure over a wide range of size scales. In 
particular, so-called "small scale structure", ranging in size from a few hundred kilometers 
up to parsec-scales, has been uncovered by a variety of observational techniques over the 
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last few decades. This structure manifests itself in neutral, ionized and molecular form 
and appears to be over-pressured with respect to the surrounding medium. Small-scale 
atomic structur e was first revealed by a VLBI HI absorption study towards the extragalactic 



source 3C 147 (iDieter et al. Ill976l ). which showed AU-scale variations in HI opacity. Other 
techniques to probe these variations include the use of optical tracers (particularly Na I D 
lines) against binary stars and globular clusters (e.g. Lauroesch et al. 1998), multi-epoch 
HI absorption measurements against pulsars (e.g. Frail et al. 1994), and mapping of HI 
optical depth against continuum sources (e.g. Faison et al. 1998). Small-scale molecular 
structure is revealed by observations of time-variability of molecular lines seen in absorption 
against quasars (e.g. Marscher et al. 1993), and small-scale ionized structures are probed 
using pulsar scintillation (reviewed by Rickett 1990). 

At present, no single theory is able to account for all forms of small-scale structure. The 
challenge is to explain the observed large ranges in pressure and temperature (e.g. Jenkins 
& Tripp 2001) in the context of a dynamic ISM. The ISM is often described by a multi- 
phase model that includes the effects of supernova explosions and X-ray ionization (McKee 
& Ostriker 1977). A feature of this paradigm is that the various phases are in approximate 
pressure equilibrium. However, it is now thought that turbulence has a significant role 
in creating a more complex medium at small scales, featuring large contrasts in density, 
pressure and temperature (e.g. Kritsuk & Norman 2002; Koyama & Inutsuka 2002; Audit 
& Hennebelle 2005; Gazol et al. 2005; Heitsch et al. 2005; Mac Low et al. 2005). A major 
obstacle is that the fundamental nature of small-scale structure is not understood. The 
discovery of "tiny" , low column density clouds (e.g. Braun & Kanekar 2005; Stanimirovic & 
Heiles 2005) suggests it could be in the form of discrete objects, but it is also argued that it 
is part of a power spectrum resulting from turbulent processes (e.g. Deshpande 2000). The 
pervasiveness of small-scale structure is also unknown. For example, a multi-epoch pulsar 
study by Frail et al. (1994) suggested that ~ 15% of the cold neutral medium (CNM) is in the 
form of AU-scale structure. However, more recent pulsar work (e.g. Johnston et al. 2003; 
Stanimirovic et al. 2003) shows structure to be more sporadic and argues that previous 
"detections" wer e likely the result o f systematic calibration errors. VLBA observations 



towards 3C 138 (iBrogan et al. I l2005l ) suggest the filling factor of atomic structure in the 
CNM is < 1 %. Thus, it is especially complicated to assess the degree to which small-scale 
structure is universal in the ISM since it is difficult to detect. 

To make progress in understanding structure generation it is important to assess the 
relative roles of dyamical processes and those involving thermal conduction and evaporation. 
For example, Slavin (2007) has argued that evaporation is unimportant in dynamically active 
regions but that it may become dominant in quiescent regions. In this paper we examine the 
stability properties of evaporation and condensation fronts separating two media of different 
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temperatures and densities. The thermal-type corrugational instabihty results from pertur- 
bations that are periodic on a surface of discontinuity, causing it to wrinkle. Its character 
has been discuss ed in the contex t of shock waves (D'yakov 1954; Landau & Lifshitz 1987) 



and combustion (IWilliams Ill964j ). The instability was first applied to combustion in Type 
la supernovae by Bychkov & Liberman (1995). In a study of the structure and stability 
of CNM/WNM phase transition layers Inoue et al. (2006; hereafter IIK06) showed that 
evaporation fronts are hydrodynamically unstable to corrugational deformations, whereas 
condensation fronts are hydrodynamically stable. It was also suggested that such perturba- 
tions could generate turbulence in the nonlinear regime w ithin a dynamical time of the ISM. 



Three-dimensional hydrodynamic numerical simulations (iKritsuk &: Norman II2002I ) suggest 
that the formation of filaments and voids in a multiphase ISM induce turbulence that can 
saturate the nonlinear thermal instability. Our primary interest here is in the potential of 
the corruguational instability to generate structure. 

The Alfven velocity in the neutral ISM is at least 10 - 100 times larger than the flow 
speed of fronts. Therefore, magnetic fields will significantly affect their dynamics. In this 
work we discuss the modifications to the hydrodynamic analysis for the scenario of a magnetic 
field orthogonal to a front, and concentrate on the sub-Alfvenic regime as the astrophysically 
relevant case. The linear effects of magnetic fields on combustion fronts in the context of 
supernovae were considered by Dursi (2004). In that work it was shown that super- Alfvenic 
flames are not significantly affected by a magnetic field, whereas in the sub-Alfvenic case 
the instability can be greatly suppressed. Dursi's work is equivalent to our evaporation front 
setup, so in this paper we extend the analysis to condensation fronts. 

We discuss the physical picture of fronts in §2 and present the linear stability analysis in 
§3. Mode analyses and dispersion relations for the super-Alfvenic and sub-Alfvenic regimes 
are presented in §4. The physical mechanism of the instability is discussed in more detail in 
§5. We apply our findings to the neutral ISM in §6, and summarize our results and mention 
future work in §7. In the appendices we present the derivation of the dispersion relation 
for the magnetized entropy- vortex wave and compare our method with that of Dursi (2004), 
to allow ease of comparison between the supernova and ISM applications. This paper is 
the first in a series exploring the linear and nonlinear effects of magnetic field strength and 
orientation on the corrugational instability. 



2. Front Structure 



We consider the physical setup of a phase transition layer (or front) separating two 
uniform media of different densities and temperatures in a 2-D plane-parallel geometry. 
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These phases result from the balance of radia tive heating and cooling due to cosmic rays and 



external radiation fields (IWolfire et al. II2003I ). The nature of a front depends on the pressure. 
There exists a "saturation pressure" such that heating and cooling are balanced inside a static 
front (Zel'dovich & Pikel'ner 1969; Penston & Brown 1970). An evaporation front exists if 
the saturation pressure exceeds the external pressure such that there is net heating across the 
transition layer. In a condensation front the saturation pressure is exceeded by the external 
pressure such that there is net cooling across the layer. The physical setup is illustrated in 
Figured! We assume the front is located in the plane x = with the ambient flow parallel to 
the X-axis and select the reference frame so that it moves with the front. A uniform magnetic 
field is orthogonal to the front, but has no effect on its structure since it does not exert any 
force and does not modify thermal conduction in the direction of the temperature gradient. 
We assume the perturbations on each side of the front are incompressible, since the ambient 
flow speed is subsonic (IIK06). 

The thickness of a front (the Field length, Xp) is given by 



wher e k(T) is the temperature- dependent thermal conductivity and n'^C is the net cooling 



rate (IBegelman fc McKee Ill990l ). For example, has the values 0.01 pc and 0.1 pc for the 
CNM and WNM, respectively (IIK06). In order to include all relevant physical processes 
one must consider perturbations with wavelengths both larger and smaller than Xp. For 
wavelengths larger than the front thickness the transition layer is treated as discontinuous 
and perturbations of the fluid dynamical variables are important. This allows for a straight- 
forward analysis but neglects the effects of thermal conduction and thermal stability (IIK06). 
To remedy this one must also consider corrugations smaller than the front thickness for which 
perturbations of the thermodynamic variables are important. In this work we only discuss 
the long wavelength, discontinous front approximation. The roles of heating, cooling, and 
thermal conduction in front stability will be the subject of future study. 



3. Phase Transition Layer Stability Analysis 

We perform a linear stability analysis for fronts undergoing corrugational deformation 
when subject to an orthogonal magnetic field. We choose a similar method to that of IIK06, 
rather than that used by Dursi (2004). We provide comments throughout this paper on the 
commonalities of all three approaches, with a more detailed comparison provided in Appendix 
B. It should be noted that Dursi 's combustion front is equivalent to our evaporation front 
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setup, in the sense that the downstream density is lower than the upstream density. Our 
work is thus an extension of his in providing the stabihty analysis for a condensation front. 

We begin by linearizing the equations of continuity, momentum conservation, and in- 



duction in the incompressible approximation to obtain 

V • (5v) = 0, (2) 

^(Pc/,d5v) + pc/,d(vu,d ■ V)5v + VP = ^ '^^) ^ ^ (3) 
Ot 47r 

^(SB) = V X [(5v X B) + (v X 6B)], (4) 



where subscripts "U" and "D" refer to upstream and downstream, respectively, p, v, P, 
and B represent density, velocity, pressure, and magnetic field, and perturbed quantities are 
denoted by the use of "5". We assume ideal MHD and neglect diffusivity. We take the 
front to be in the plane x = and assume perturbations in the (k, B) plane of the form 
exp[i{kxX + kyU — cut)], so that unstable modes have Im(ti;) > 0. Substituting the form of 
the perturbations into equations ([2]), ([3]) and (jl]), and taking B = i?x (i.e. orthogonal to the 
front) yields the following equations: 

kjv^ + ky6Vy = 0, (5) 

pu,D{-uj6vr, + kxVu,D5vx) = -kxSP, (6) 
Pu,D{-^^Vy + kxVu,D^Vy) = -ky6P + — {kJBy " kySB^), (7) 

- ui6Br, = ky{vu,DSBy - B6vy), (8) 

- uSBy = k,j,{B5vy - vu,DSBy). (9) 

The perturbations are decomposed into two surface waves that are localized to the 
front, and two magnetized entropy- vortex waves that are advected by the flow. We adopt a 
notation similar to that of IIK06 and denote the surface waves upstream and downstream of 
the front by 6vl ^ and Svy^\ with dispersion relations given by A;i ■* = —iky and k^^ = iky, 
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respectively. The magnetized entropy- vortex waves are represented by Sv'^ ^ and 5v\^^\ 
Their dispersion relations are given by kx~'^ = ujjiv — va) and k^x^^ = uj/{v + va), where 
the Alfven speed is given by va = B/y/Anp. The nature of these magnetized entropy- vortex 
waves and the derivation of their dispersion relations are discussed further in Appendix A. 

The directions of propagation of the magnetized entropy- vortex waves depend on whether 
the flow is super- or sub-Alfvenic. In the super-Alfvenic case no magnetized entropy-vortex 
waves can propagate upstream, so Svlf ■* and Svy^~^^ are both downstream perturbations. In 
the sub-Alfvenic regime 6vif~^ propagates upstream whereas Svy^~^^ is a downstream pertur- 
bation. In §4 we derive separate characteristic equations for each regime. 

The X— component velocity, magnetic field, and pressure perturbations can be rewritten 
in terms of the component velocity perturbations upstream and downstream of the front, 
using equations ([5]) - ([9]). We represent perturbed quantities as sums of the surface waves 
and magnetized entropy-vortex waves: 



H = E''4°' (10) 



*'^ = -ET|yH°'. (11) 

(12) 

SBx = y 6it\ (13) 

a [Vu,Dkx - Uj) 



B'iki-'^' + kl) p^,^ 



{vu,Dk^^ - Uj) 



(14) 



A'Kky{vu,Dk^x^ - u) K 

where a corresponds to the different types of perturbation: (-),(-!-), (s-), and (s-|-). These 
quantities are computed upstream and downstream of the front, with the appropriate per- 
turbations substituted for 6v'^\ depending on whether the flow is super-Alfvenic or sub- 
Alfvenic. 

To proceed, we require the jump conditions across the front. The normal and tangential 
unit vectors are given by: 

d . d 

n = x-y—6xf, t = y + x—Sxf, (15) 
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where 5xf is the front displacement. These are used to compute the velocity perturbations 
normal and tangential to the front: 

5vn = Svx + iuj6xf, 6vt = Svy + ikyVu^D^xj, (16) 

(IIK06). The expression for 6vn includes the motion of the front itself. 

The jump condition due to the continuity equation is given by 

[pSun] = 0, (17) 

where represents the jump in quantity \E' across the front. Following IIK06, we modify 
this condition by assuming that the structure of the front is unchanged by a perturbation 
such that: 

PuSVn,U = PDSVn,D = 0. (18) 

This is equivalent to the condition used by Dursi (2004) that both sides of the front travel 
with the same velocity {[6vx] = in our notation). 

The V ■ B = condition implies that the normal component of the magnetic field 
perturbation is conserved across the front, thus: 

[SBx] = (19) 

Linearization of the normal component of the momentum flux conservation equation yields 

2p6vr.v^ + 6P- = 0. (20) 

4:71 

By invoking the continuity of 6Bn the corresponding jump condition is given by: 

2pv5vn + 5P\ = 0. (21) 

Linearization of the tangential component of the momentum flux conservation equation yields 
the jump condition 



pvovt ■ 



4. »• 

where 6Bt = ikyB6xf + 6By. By substituting the expressions for the velocity, pressure and 
magnetic field perturbations, equations (fTTj) - (fT^ . into the jump conditions, equations (ITHl) 
- (122|) . we obtain the characteristic equation of the system. This is written in the form: 



Mt] = 0, (23) 
where the matrix M is derived from the jump conditions and 

V = {6xf, 6vl~\ 5v^^\ 6vl;-\ 5^^+))^, (24) 

To obtain the dispersion relation governing the perturbations we set the determinant of M 
equal to zero and solve for uj. 
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4. Mode Analysis 



In this section we provide the characteristic equations for super-Alfvenic and sub- 
Alfvenic fronts, in terms of the density jump across the front, = po/ Pui and the ratio of 
the downstream Alfven speed to the downstream flow speed, A = vad/vd- Super-Alfvenic 
fronts have Vu > vau and vd > vad, i-e. \/r^ > A on the upstream side and 1 > A on the 
downstream side. Sub-Alfvenic fronts have vu < vau and vd < vad, i-e. \/r^ < A on the 
upstream sid e and 1 < A on the dowstream side. A third regime, called "trans-Alfvenic" 
( iDursi II2OO4I ). may be identified in which the flow speed exceeds the Alfven speed on one 
side, but is less than the Alfven speed on the other side, i.e. A takes a value between ^/r^ and 
unity (v^^ < A < 1 for evaporation, 1 < A < ^/r^ for condensation). The front behavior in 
this regime cannot be found using the methods discussed in this work, as it is possible for 
more waves to propagate than there are boundary con ditions const raining them, and thus 
no unique solution exists for the initial- value problem (IDursi II2004I ). We argue in §6 that 
the sub-Alfvenic regime is of most relevance to the interstellar medium, so do not provide 
any further discussion of the trans-Alfvenic regime. The stability properties for both front 
types in the various regimes are summarized in Table 1. 



4.1. Super-Alfvenic Fronts 

If a front is super-Alfvenic no entropy-vortex waves can propagate upstream, so the 
total upstream y-velocity perturbation is simply given by the upstream surface wave: 

6vy = (25) 

Downstream, the total y-velocity perturbation is given by: 

Svy = 5vl+^ + 5vl;-^ + 5vl;+\ (26) 

The velocity, magnetic field, and pressure perturbations on each side of the front are written 
in terms of these total y-velocity perturbations, and substituted into the jump conditions. 
Collecting terms in each perturbation, the matrix, M, of the characteristic equation for the 
super-Alfvenic case is given by: 

/ u -1 \ 

a; 1 ^(1-A) ^(1 + A) 



M 



1 + ^ ^(l-^)(2-^) '-^{1 + A){2 + A) 

tujn(rd-l) + " -l + A " -I -A 

. \-\ 1 + i , 

(27) 
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where 00 d — kyVo- 

Taking the determinant of this matrix yields the dispersion relation: 

: — : — r — U (ze) 

ArdUJoird^D - i(^){(^D + i'^)'^ 
This can be simphfied by removing two "trivial" modes that do not cause a displacement of 
the front. These modes are given by a; = i{l ± A)u;d- Factoring these out, the dispersion 
relation can be written as: 

^3 + tVdUJn^^^y + 7^^(1 + 2A' - 3rd)u + ^r^^i,^^^ = 0. (29) 
[rd + 1) [rd + 1) (1 + rd) 

In the limit A — > the hydrodynamic dispersion relation of IIK06 is recovered. The super- 

Alfvenic stability behavior is similar to that of the hydrodynamic case: condensation fronts 

remain stable while an unstable mode is present in evaporation fronts. As found by Dursi 

(2004), the growth rate of this unstable mode is slightly enhanced. We discuss the physical 

origin of this behavior in §5. 



4.2. Sub-Alfvenic Fronts 

For sub-Alfvcnic fronts the upstream y-velocity perturbation is composed of a surface 
wave and entropy-vortex wave: 

5vy ^ 5v^-^ + Sv^;-\ (30) 
and the downstream y-velocity perturbation is given by: 

Sy^ = Sy(+) (31) 

The matrix, M, of the characteristic equation for the sub-Alfvenic case is: 

/ ^ -1 '^V^diV^d-A) \ 

u 1 ^(1 + A) 

M= 1 + ^ 1-^ ^^^(A-V^J(2-^) ^^{1 + A){2 + A) 

iujD(rd-l) l-^^i^ -i + i^^ -1-A 
— — — "'^ :/ld_i 1 + - 

(32) 

There are two trivial modes, given by a; = (1 + A)u}£i and a; = i{\/r ^A — rd)i^D^ which do 
not displace the front. Factoring these out, the dispersion relation may be written as: 

1 + rd 1 + rd 1 + r<i 

(33) 
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The stability properties are the opposite of what is obtained in the super-Alfvenic case: 
the evaporation front is stabihzed, but an unstable mode exists in condensation fronts, 
the growth rate of which decreases with increasing magnetic field strength. For very large 
magnetic field strength, as is expected for the ISM (see §6), the growth rate is approximately: 

u ^ ' ^^^^ ^ un (34) 

This can be written in more physical terms if we assume that material crosses the front in 
a time of order the cooling time, given by tf. = 5kT/2nA, where nA is the cooling rate. For 
our calculations we use the cooling function of IIK06, in which nA = 2.0 x 10~^^ ergs s~^ 
corresponds to thermal equilibrium. Assuming the most unstable scale is twice the Field 
length (as is found by IIK06), the growth rate is given by: 

_ zry'(r,-l) 27rnA 
2A^ 5k^ ^^^^ 

The growth rates of the unstable condensation mode in the sub-Alfvenic regime are plotted 
as a function of A for various in Figure O The role of the density jump across the front 
in driving the instability decreases in importance with increasing A. 

We have also examined the eigenfunctions of the sub-Alfvenic matrix to assess the 
contributions of the various perturbations. The amplitudes of the entropy-vortex waves 
relative to the surface waves are of order A~^. Thus, the instability produces very little 
disturbance upstream or downstream and remains localized to the front. 



5. Physical Mechanism 

To investigate the mechanism of the corrugational instability we begin with the simple 
case of a front in the hydrodynamic limit, in which evaporation fronts are unstable whereas 
condensation fronts are stable. Since the normal component of the velocity perturbation is 
zero on both sides of the front, we use equation (fT6ll to write: 

Svx = —iuSxf. (36) 

If iuj is a positive, real number, as is the case for evaporation fronts, the x-velocity pertur- 
bation will reinforce the displacement of the front and continue to drive the instability. In 
condensation fronts 6vx opposes the displacement and thus stabilizes the front. Since all 
growing modes have positive imaginary u, this argument is true in general for all unstable 
fronts, whether super-Alfvenic or sub-Alfvenic. 



- 11 - 



Streamline curvature can also be used as a diagnostic of stability, as has been done for 



flame fronts in Type la supernovae (jPursi II2004I ). In that work it was argued that in regions 



where the streamlines "fan out", the local flow velocity into the front is reduced and the 
front can propagate further ahead. We demonstrate this analytically, by calculating the angle 
formed between a given streamline and the front normal, to first order in the perturbation 
amplitude, using: 

nx vu,D Svt 
sm eu,D = n = TT (37) 

\v\ \v\ 

In the hydrodynamic case, Svt is continuous across the front, so we have sin 9u = Svt/vu 
and sin 60 = Svt/vo. Evaporation fronts have vr, > vu, so the downstream angle is smaller 
than the upstream angle so the streamlines bend towards the front normal. In the case 
of a condensation front vu > vd so the streamlines bend away from the normal. This is 
illustrated in Figure [31 In regions where the streamlines diverge the x- velocity perturbation 
is negative {Sv^ < 0) and so the amplitude of the corrugation is enhanced. 

In the case of a magnetic field orthogonal to the front, this analysis is more complicated, 
since the condition on 6vt is changed. We derive it using the y-component of the linearized 
induction equation f|T2|) and the tangential component of the linearized momentum equation 
f l22]) . We assume a; << k^vu^Di which can be justified a posteriori, and write the tangential 
momentum balance condition as: 

- 5vtu = (1 - A^) SvtD. (38) 

Upstream and downstream angles may be compared by writing equation (l38l) as: 

6vtu — 1 1 5vtD , N 

(39) 



vv — - Ird vd ' 

rd 

such that 

- 1 

sin Ou = sin Od-t^ • (40) 

In the super-Alfvenic regime we consider A to be small compared to both 1 and r^, so 
equation (HO!) becomes sinOu ^ r^^ sin6'£)(l — A^ + A^/r^). If < 1, the change in streamline 
angle across the front is increased by the magnetic field, so the instability is enhanced in 
evaporation fronts, as noted in §4.1.. For > 1 the angle decreases across the front so 
condensation fronts are stable. In the sub-Alfvenic regime, in which A » 1 we obtain 
sin6u ~ sin6'£i[l + (r^ — l)/^^], which is consistent with the stabilization of evaporation 
fronts and the slight destabilization of condensation fronts. Thus, the streamline picture 
explains the regimes of stability and instability found from the dispersion relations given by 
equations fl29l) and fl33l) . 
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6. Application to the ISM 



The structure of fronts is calculated in 11K06 by solving the energy equation while 
assuming thermal equilibrium at infinity. In future work we will study the effects of magnetic 
fields on front structure, but for now we simply use their results to examine the growth rates 
of unstable modes for parameters appropriate to the neutral ISM. The magnetic field strength 
parameter, A, can be written in terms of the upstream gas conditions as: 



where R is the molar gas constant, fi is the mean molecular weight (assumed to be unity), Tu 
is the temperature in the upstream gas, and Pu is the upstream pressure. The parameters 
corresponding to examples of steady state evaporation and condensation fronts, derived from 
the front calculations of IIK06, are presented in Table 2. Both types of front are sub-Alfvenic, 
with their dynamics dominated by the magnetic field. In this regime evaporation fronts 
are stable, so we now calculate the growth rate of the unstable condensation mode. Using 
equation (l35l) and the parameters given in Table 2, we find a growth timescale tgrow = l/cu ~ 
1.8 X 10^^ s. This would imply that the corrugational instability is essentially stabilized in 
the neutral ISM. 

However, if condensation proceeds more rapidly the instability may grow within a dy- 
namical time of the ISM. Inoue et al. (2007) argue that the condensation speed during the 
growth of a thermal instability that can generate clouds is of order 1 km s~^. If we use 
this velocity with the other parameters in Table 2 we obtain A ~ 360. Assuming the most 
unstable scale is of order 0.1 pc we find a growth timescale tgrow ~ 0.6 Myr. This would 
imply that the corrugational instability has a significant effect on cloud dynamics. 

Since we have neglected the thickness of the front in this analysis we cannot provide fur- 
ther comments on the physical scale of the instability or its growth time. The full significance 
of the unstable condensation mode and the actual speed of the condensation process will only 
become apparent through completing the finite front analysis and nonlinear simulations. 



We have shown that for the case of a magnetic field orthogonal to a phase transition 
layer unstable modes can exist in both evaporation and condensation fronts, depending on 
whether the flow is super-Alfvenic or sub-Alfvenic. This is in contrast to the hydrodynamic 
case in which an unstable mode only occurs in evaporation fronts. We demonstrated this by 



A 



VAD B RTuTd 




(41) 



7. Summary and Conclusions 
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deriving linear dispersion relations, given by equations fl29l) and fl33|) . and from an analysis 
of streamline geometry near the corrugated front. 

In the super-Alfvenic case the instability in evaporation fronts is enhanced from the 
hydrodynamic case, which is apparent from the increase in streamline bending. In the 
sub-Alfvenic case evaporation fronts are stabilized and it is the condensation front that is 
unstable. A simple formula for the maximum growth rate of the condensation front instability 
is given by equation (^S^, from which it is shown that under conditions appropriate to steady 
fronts in the neutral ISM the growth time is sufficiently long that such systems are effectively 
stable. We make the simple argument that the corrugational instability is suppressed by the 
rigidity of the magnetic field. However, if condensation proceeds sufficiently rapidly the 
growth timescale may become less than the dynamical time, such that the instability is 
important for cloud formation dynamics. 

Forthcoming papers in this series will consider perturbations smaller than the front 
thickness for both orthogonal and tangential magnetic field orientations. The full signifi- 
cance of the unstable modes will only be revealed through studies of their nonlinear behav- 
ior. Numerical simulations will be performed to assess the nonlinear characteristics of the 
instability and its potential for generating small-scale structure in the ISM. 

We acknowledge support from NASA ATP Grant NNG05IG09G, and NSF Grants 
AST-0507367 and PHY-0215581. We are very grateful to our anonymous referees whose 
recommendations greatly improved this paper. This work has also benefited from useful 
discussions with F. Heitsch, K. M. Hess, A. S. Hill, S.-I. Inutsuka, N. A. Murphy and S. 
Stanimirovic. 
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APPENDIX A 



Derivation of the Magnetized Entropy- Vortex Wave Dispersion Relation 

The magnetized entropy-vortex wave is a vortical perturbation that bends magnetic 
field lines and generates an Alfvenic response. The derivation of its dispersion relation is 
provided here for the interested reader. We seek a perturbation downstream of the interface 
satisfying V ■ (5f = while including t he effects of magnetic field s. The derivation is modeled 



after that for the unmagnetized case (ILandau fc Lifshitz 1119871 ). 



We begin by writing linearized equations for the entropy, pressure and magnetic field 
perturbations: 

i{kxv - Lj)5S = 0, (A-1) 



i{kxv - uj)5P = (A-2) 

i{kxV — uj)5Bx = —ikyBSvy, (A-3) 

i{kxV — uj)6By = ikxB6vy. (A-4) 

The velocity perturbations may be found using the x— and y— components of the 
linearized momentum equation: 

SP 

i{kxV — ol!)5vx = —ikx — , (A-5) 

P 

SP iB 

i{kxV — uj)5vy = —iky h - — {kxSBy — kySBx). (A-6) 



The energy and linearized induction equations (A-2, 3, 4) can be used to rewrite these 
components as: 

Sy^ = k,cl^J^4?^±tytyl^ (A-7) 



= +ky5Vy) ^ _ 
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Multipying flA-7|l by and flA-81) by ky and rearranging yields: 

k'^cl 1 . , . . k'^v\ 



1 



{k^v - 



{kJVx + kydVy) = kySVy. 



(A-9) 



{k^v — uj 

The only solution that gives V-5v = is /c^ = 0, i.e. kx = iiiky, corresponding to exponen- 
tially decaying perturbations. These are surface waves. 



We proceed with the derivation by taking the linearized equations of motion and induc- 



tion: 



p 



6P + 



B6B 

An 



B V6B 



d_ 

di 



5B = Vx(5v X B), 



and computing the curl of these to obtain: 



d_ 

di 



|(Vx^v) = ^(Vx^B). 



(Vx5B) = VV-(5v X B) - V^((5v x B) 



The 2— component of the curl of the induction equation is: 
-2cu(V X 5B), = V^'°^- ' - '^^ 

and the vorticity is: 
Combining these yields: 



- 2cu(V X 6B), = V\B6vy) = -k^B6vy, 
(V X 5v). = ^6vy - ^6v. = '^Jvy. 



'ik"^ 
iuj{ - — 5vy 



zk^B f-ik^B , 

Of. 



Arrp 



UJ 



(A-10) 
(A-11) 

(A-12) 

(A-13) 

(A-14) 
(A-15) 

(A-16) 



The solutions to the above are fc^ = 0, or uj"^ = k^v\, so in a moving medium the 
dispersion relation of the magnetized entropy-vortex wave is given by: 



{u - k^v) = k^v^ 



where 



kx 



(A-17) 



(A-18) 



This can be compared to the dispersion relation for the unmagnetized entropy-vortex wave, 
which is simply kx = -■ 
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Appendix B 



Comparison with Corrugational Instability in Burning Fronts (Dursi 2004) 

The effect of a perpendicular magnetic field on the linear instability of flame fronts 
in Type la supernovae was demonstrated by Dursi (2004). Here, we compare his method 
with ours and summarize his findings, to allow comparison between the ISM and supernova 
applications. The physical setup of the supernova flame front is similar to that shown by the 
left side diagram in Figure [H one should consider "Phase I" to be unburned fuel and "Phase 
II" to be burned ash. Dursi's x-z coordinate system is obtained by rotating ours 90 degrees 
about the y-axis and then 90 degrees about the z-axis. The magnetic field is orthogonal to 
the flame. 

Dursi assumes incompressibility and that the reaction driving a supernova flame front 
is sufficiently exothermic such that the density of the burned material is less than that of the 
incoming fuel. Thus, his scenario is equivalent to our evaporation front setup. To allow ease 
of comparison we rewrite his notation in terms of ours. The density ratio of the unburned 
to burned material is given by: 

a = ^ = - (A-19) 

Pb Td 

The strength of the magnetic field is represented by the ratio of the Alfven speed in the 
unburned material to the fiow speed: 

a„ = ^ = ^ (A-20) 

The growth rate itself is given by: 

n —iuo 

n=-— = A-21 



It can be shown that the jump conditions we have derived and those employed by Dursi 
are the same. The only differences are our assumption that the front structure is unchanged 
by perturbations and Dursi's additional condition that each side of a front travels at the 
same velocity {[5v,^ = 0). Although this is not a condition that we directly impose it is 
still satisfied by our eigenfunctions. The dispersion relations we present are identical to 
those that are derived in Dursi (2004) and can be recovered by substitution of the relations 
given above. We have extended Dursi's work by giving full consideration of the behavior of 
condensation fronts and applying our results to the neutral ISM. 
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upstream downstream 
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Fig. 1. — Physical setup of ISM phase transition layers undergoing corrugational deformation in a 
plane-parallel geometry. The regions to the left and right of fronts are identified as upstream and 
downstream, respectively. The phases are uniform and have different temperatures and densities. 
Vu and are the upstream and downstream flow velocities, respectively. In an evaporation front 
there is net heating across the transition layer and in a condensation front there is net cooling. A 
magnetic field is orthogonal to the front. 
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12 3 4 

log A 

Fig. 2. — Growth rate of the unstable mode in condensation fronts of various as a function of A 
for the sub-Alfvcnic regime. The growth rate is given in units oi ujd. As the magnetic field strength 
is increased, the role of the density jump across the front in driving the instability becomes much 
less significant. For A > 1000, as expected in the neutral ISM, the growth rate is very slow, owing 
to the rigidity of the magnetic field. The blacked-out area in the top left corner covers the region 
of parameter space {A < ^/r^) that is outside of the sub-Alfvenic regime. 



UNSTABLE FRONT 



STABLE FRONT 




Fig. 3. — Upstream and downstream streamline angles for stable and unstable fronts. The front 
is represented by the thickest line and the dotted line is normal to it. Ou and 6^ represent the 
upstream and downstream angles, respectively. The thin solid line represents a given streamline. 
In an unstable front the downstream angle is smaller than the upstream angle, so streamlines are 
bent towards the normal. This accelerates the front, enabling it to propagate further ahead and 
drive the instability. In a stable front the streamlines diverge and the corrugation is smoothed out. 
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Table 1. Summary of Front Stability Properties 



Front Type 


Hydrodynamic 


Super-Alfvenic 


Trans- Alfvenic'^ 


Sub-Alfvenic 


Evaporation 


unstable 


unstable 


unknown 


stable 


Condensation 


stable 


stable 


unknown 


unstable 



^It is beyond the scope of this work to ascertain front stability in this regime. Non- 
linear simulations will be required. 



-23- 



Table 2. Example parameters for steady state evaporation and condensation fronts in the 
ISM (calculated from IIK06), assuming B — 6 /i G 



Front Type 


P (K 




nu (cm ^) 


rd 


vu (km s ^) 


A 


Evaporation 


1200 


35.5 


33.8 


0.004 


0.14 


0.95 


C ondensat ion 


4000 


8115 


0.49 


371 


0.0302 


11,900 



